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We consider the generation and evolution of quantum fluctuations of a massive 
nonsymmetric gravitational field (-B-field) from inflationary epoch to matter era 
in the simplest variant of the nonsymmetric theory of gravitation (NGT), which 
consists of a gauge kinetic term and a mass term. We observe that quite generically 
a nonsymmetric metric held with mass, tub — 0.03(////10^'^ GeV)^ eV, is a good 
dark matter candidate, where Hi denotes the inflationary scale. The most prominent 
feature of this dark matter is a peak in power at a comoving momentum, k ~ 
^/mBI^o/{l + ^^cq)"*^^^, where Zgq is the redshift at equality. This scale corresponds 
roughly to the Earth-Sun distance. 

PACS numbers: 



I. INTRODUCTION 

Einstein's theory of relativity has passed all experimental tests on laboratory and inter- 
mediate scales. On cosmological scales however, in order to get a successful description of 
the Universe's dynamics one requires addition of dark matter and dark energy, both of un- 
known composition. An alternative is to extend Einstein's theory and hope that a modified 
theory of gravitation would describe gravitation on cosmological scale without a need for 
dark matter and/or dark energy. Examples of extended theories of gravitation which have 
both been used as alternatives to dark matter are the nonsymmetric theory of gravitation 
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(NGT) 1^ and MOND 

In Ref. jsj we assumed that the antiymmetric tensor field B is dynamical, and consid- 
ered the cosmological evolution of quantum fiuctuations generated in de Sitter inflation of 
a massless 5-field, which is related by a duality transformation to the Kalb-Ramond axion 
field. We have further shown that the evolution of a massless B field mimics that of a 
very light field. We then presented some early results on the evolution of a massive 5-field, 
and observed that the scaling changes when the field becomes nonrelativistic in radiation 
era. From our analysis of quantum fiuctuations of the physical mode of the antisymmetric 
field, which corresponds to the longitudinal 'magnetic' component, it follows that the field 
couples conformally during de Sitter infiation, such that the i?-field correlator exhibits the 
spectrum of conformal vacuum fiuctuations. This is contrary to what has been claimed 
based on studies of the Kalb-Ramond axion, which couples to gravitation just like a mass- 
less minimally coupled scalar field, and therefore exhibits a nearly scale invariant spectrum 
during infiation. In fact, because of the conformal coupling during infiation, there cannot 
be a physical observable which exhibits a scale invariant spectrum during infiation or sub- 
sequent epochs. There is amplification however, which is induced by the matching at the 
infiation-radiation transition, and based on which the spectral amplitude gets enhanced at 
superhorizon scales with respect to the conformal vacuum, but not enough to get a scale in- 
variant observable. Finally we observed that the spectral features of the B field fluctuations 
produced in inflation are similar to that of gravitational waves, making it thus an alternative 
probe of inflationary scale. There is an important difference in the spectrum between the 
primordial gravitons and the B fleld fluctuations: while the primordial gravitons exhibit a 
scale invariant spectrum on superhubble scales, the -B-fleld fluctuations are suppressed. 

Building on an earlier work of Damour, Deser and McCarthy and Clayton [7], Janssen 
and Prokopec ^ have recently shown that there is no nonsymmetric geometric theory which 
yields a dynamical B field and which is "problem free," in the sense that the 5-field dynamics 
is ghost free and it is not marred by instabilities. Furthermore the authors of js| have shown 
that the most general (problem free) quadratic action in the B field is of the form, 



S = She + Sb, Sb = J d xCb 
C-B = -V^g^H^^.H^'^P - v/^i (m| + in)B,,B^'' , (1) 
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where 

H^iyp = d^B^p + dyBpp + dpBpy (2) 
is the antisymmetric field strength and 



HE 



/A£„.. £„, = -^(k + 2a) (3) 

is the Hilbert-Einstein action, with A being the cosmo logical constant. In this paper we 
work in the units in which, l/(167rG'Ar) = 1. From Eq. we see that the coupling to the 
Ricci scalar TZ is the only allowed coupling to curvature. 

In order to simplify the problem further, we drop the coupling to the Ricci scalar and set 
^ = in Eq. (^. This term is expected to act as an additional mass term during infiation 
era and thus introduce additional breaking of conformal invariance. During radiation and 
matter era however, where the Ricci scalar is either zero (radiation era) or of the order of 
the Hubble parameter (matter era), we expect that the Ricci scalar term has no significant 
impact on the evolution of the i?-field. 

In a geometric theory the mass term is naturally induced by the cosmological constant, 
in which case j8|, 

m| ~2A(l-2p + 8(T), (4) 

where p and a are defined by the following decomposition of the metric tensor, 

9,u = g^u + B^, + pB^B + aB^gp, + 0{B%) . (5) 

This is a general metric tensor decomposition in terms of the symmetric metric tensor {gpv) 
and antisymmetric metric tensor {Bpy). Yet there is no reason to assume that the mass 
term is fully of geometric origin, and thus the relation (j3)) needs not to hold in general. 

Equation (0)) is also significant because it implies that, if the B field is of a geometric 
origin, then it would be unnatural to assume that its mass vanishes. Indeed the current ob- 
se^fo... based on the ..ni.osity.edshift .elation of distant supe.novae la flQ. suggest 
that the cosmological term today is of the order, A ~ 10~*^ GeV^. To maintain generality 
we assume in this work that the mass is unspecified and study its cosmological impli- 
cations. In order to do that, we canonically quantise the massive B field in infiation, and 
follow its subsequent dynamics in radiation and matter eras. Our main finding is that the 
B field with a mass of the order, 

mg-0.03( ^ j eV (6) 
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which corresponds to a lengthscale, m~j^ ~ 7 x 10~^(i///10^^ GeV) m, is a good dark 
matter candidate. 

Since the B field is produced in infiation and does not couple to the matter fields, its 
spectrum is highly nonthermal. Indeed, we find that the spectral power is peaked at a 
comoving momentum, k ~ v^mB-f^o/(l + -2cq)^^^) where 2;eq is the redshift at matter-radiation 
equality, and corresponds to a physical scale at structure formation {z ~ 10), ^phys ~ 
2 X 10'' km. This peak is generated as a consequence of a different nature of the vacuum 
states in inflation and radiation era. This is the main feature by which this dark matter 
can be distinguished from other dark matter candidates, which typically obey a thermal 
statistic. 

Another important feature are Fourier space pressure oscillations, which occur after the 
second Hubble crossing. Although we flnd that the pressure of the B-fleld drops to zero 
before the decoupling of the cosmic photon fluid, the Fourier pressure components exhibit 
signiflcant oscillations. These oscillations may have a potentially observable impact on the 
gravitational potentials, and they are thus the second distinct feature of our dark matter 
candidate. The physical signiflcance of these spectral pressure oscillations should be further 
investigated. 



II. STABLE LINEARISED ACTION 



Taking the nonsymmetric gravity theory (NGT) as a starting point, one can perform an 
expansion of some general covariant geometrical action, in terms of antisymmetric pertur- 
bations on a symmetric background metric. The (linearised) action then reads 

^NGT = j d^x^ {R-2A- ^F.^pF^'^P 

+9:,B^pB''m^p + 0{B^)], (7) 



where the parameters 61,62, 63 and tub are deflned by the initial action. We choose to work 
up to second order in the nonsymmetric theory, and we will raise and lower indices with the 
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symmetric background FLRW metric, 

g,ug'' =K. (8) 

g^.B'^'^ =B/. (9) 



In Ref. p], it was argued that NGT in general contains propagating ghosts or unacceptable 
constraints on dynamical degrees of freedom, when no cosmological term is present. For a 
simple choice of the action, 6'„ = = {n = 1,2,3), no ghosts are found in curved 
backgrounds j^. Nevertheless, when considered in a FLRW or Schwarzschild background, 
the theory with O-l^ and/or 7^ can develop instabilities even in the presence of a 
cosmological term [q], putting into question the geometric origin of the antisymmetric tensor 
field. In the light of these results, we choose 6i = 62 = 0. In addition for simplicity we choose 
63 = 0. The linearised action considered in this paper hence reads 

Sue + SB=j d^x^ {R-2A- ^.F.^pF'^'P + im|5^,5^'^} . (10) 

In this equation we recognise the antisymmetric field strength of an antisymmetric 2-form, 
analogous to the Kalb-Ramond field |lll |. 

FfjLPp 9^Bi,p -\- diyBpp -\- dpB^u 

= V^B.p + V,Bp^ + VpB^,. (11) 



A. Field equations 

Following the action (fTUj) . the field equations for B^^ read 

a{r]fd''Fpp, - ^Fo^, + mla{r]fBp, =0. (12) 

If we take the divergence of equation (fT^. multiplied with a factor of airf)^'^, we are led to 
a self-consistency condition, 

mld^B^, =0. (13) 

Analogous to Maxwell theory, the antisymmetric tensor Bp^ can be decomposed in electric 
and magnetic components, 

Boi = —BiQ = Ei (14) 
Bij = —Bji = —eijkBk- (15) 
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In this decomposition the field equations become 

E^{x) =0, (16) 
{d^ + mlairiY)E^{x) =0, (17) 

(d' - + mlair^r^ B\x) =0, (18) 

doET - V X 5^(x) =0. (19) 

Here we use superscripts ^ and ^ to denote the transverse and longitudinal parts of the 
vectors E and B. Apparently the theory in this form possesses three degrees of freedom. 
However, if the mass term disappears, the field equations reduce to 

dPa{7^)-^Fp^, = 0. (20) 

Equation (fT!^ no longer results as a self-consistency condition. Though, the theory has 
gained a gauge freedom, 

such that equation (fT!?|l may be imposed as a choice of gauge. Equation is then by 
definition solved by 

aivY^Fp^, = ep^,,9>. (21) 

Hence, the massless theory has only one (pseudoscalar) degree of freedom, which is known 
as the Kalb-Ramond axion. One could also have made an analysis of the degrees of freedom 
of this theory by the means of its dual action, which at linearised level reduces to a massive 
Abelian gauge field action, and therefore has three physical degrees of freedom How- 
ever, in the presence of sources the dual theory contains nonlocal source terms, justifying 
the analysis in terms of the present variables. Note that the dual theory of the massive 
antisymmetric tensor field should not be mistaken for the massive Kalb-Ramond axion. 



B. Power spectrum 

What we eventually want to calculate is the spectrum of matter density fiuctuations, 
Vspik,r]), defined as p^^ {6p{x,r])6p{x' ,r])) = J ^Vsp. In this case Vsp = ^ |%f , and 5^ 
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is the fourier transform of ^. po denotes the average energy density and (■) denotes the 
averaging procedure. In the case of NGT, the power spectrum (i.e. the energy per decade 
in momentum space) is given by 

VB{k,v) = ^To'{k,v), (22) 
where T^,y is the stress-energy tensor of NGT, 

= -L-^5b. (23) 



fjLV 



That is, 



+a{r])^ml (^E{r], f )' + 5(r/, f)^) | . (24) 

III. CANONICAL QUANTISATION 

Due to the antisymmetry we find that the Oi-components of the canonical momenta are 
identically zero, 

^oi ^ Vzlpooi = 0. (25) 



This residual gauge freedom is fixed by adding the Fermi term 

C^C + ^ iVpEP^") V,5"^ (26) 

such that the canonical commutation relations can still be imposed in a covariant way. The 
Fermi term does not affect the longitudinal magnetic component of the field. 

In the next section we discuss the dynamics of the field during inflation, radiation and 
matter era. One result is that all momenta of interest, during inflation correspond to physi- 
cal momenta which are highly relativistic. Effectively the theory during inflation is massless, 
hence the transverse degrees of freedom cannot obtain any energy. From the classical field 
equations for the transverse components, there is no indication that these degrees of free- 
dom should become important at any time. In the following we will therefore neglect the 
transverse components. 
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In the context of quantising the theory with the Fermi term, it is then sufficient to say 
that the theory can be covariantly canonically quantised, and that for the component of 
interest we have the commutation relation 



UB,{x),BL{y) =--8\x-y). 



(27) 



Note that in Ref. p| it is shown that in fact the transverse degrees of freedom contain 
the instabihties of the geometrical theory, involving couplings to the curvature tensors in 
the action ((71). Canonical quantisation of transverse degrees of freedom is carried out in 
Ref. Q- 



IV. INFLATION, RADIATION AND MATTER ERA 



We can perform a Fourier transformation 

B'^(x) 



with the equal time relations 



=0, 
=0, 



The field equation for in momentum space becomes 

dl + + mlaiv)'^ B\v, k) =0, 

Equations are satisfied if the Wronskian satisfies 



W[B'^{ri,k),B*^{ri,k)] =ia{ri) 



We can conformally rescale B^ 



B' 



a{r]) 



B' 



(28) 
(29) 

(30) 
(31) 
(32) 

(33) 

(34) 

(35) 
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such that we have the field equation for 



a"{ri)a{ri) -2a'{riy 
a{r]y 



dl + k^ + m\a{rif + 



=0. 



(36) 



The expression for the power spectrum ()22|) containing only B^ in momentum space becomes, 



a' ^ 



doB'^ik,r]) + -B'^{k,r]) 
a 



+ {k' + a'mi,) B'^ik.T]) 



(37) 



A. De Sitter inflation 



Consider De Sitter infiation, in which the scale factor as a function of conformal time r) 
is given by 



a{r]) 



-1 



(with — oo < r] < ^' 



with Hj the Hubble constant during infiation. The field equation reads 



solved by 



Here zi""^ denotes any general pair of Bessel functions |13[ that forms a basis for the solution 
space. For every momentum k there has been an t] during infiation such that fcphys > mBa{ri). 
As during infiation the effectively massless B^ field is conformally invariant, its vacuum state 
is given by the conformal vacuum jl^ . at early time (as r] — oo). The vacuum state 
mode function reads, 

'-^V rr{2) 



di + e + 



B 



B^^fiv,k) =0, 



,(2) 

2 " 



{kr]). 



(38) 



(39) 



(40) 



Brnf{v,k) 




4 ^i^l-(4m|///|) 



(kri). 



(41) 



Note that in the limit when ^ 0, the field is in conformal vacuum. This is the relevant 
limit, since we are interested in the physical momenta which at the end of infiation satisfy 
the relation, Hj ^ k/a{riend) ^ ^b- Indeed, the physical momenta that are within the 
Hubble radius today (and of course smaller than the comoving Hi) all satisfy this condition. 
Moreover, since conformal invariance is broken by the stress energy tensor ()24j) . from the 
point of view of energy density, superhorizon modes in fact undergo amplification during 
infiation. This is in agreement with the stress energy tensor of the massless Kalb-Ramond 
axion, which to a good approximation corresponds to the dual description during infiation. 
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B. Radiation era 

We choose the simple model of a sudden transition from inflation era to radiation era, 
with a scale factor 



a{r]) = HjT] with < r] < 7]^). 



(42) 



Conformal time at the time of radiation- matter equality is denoted by rje- 
The field equation for the longitudinal field in this era is given by 



+ e- -do + mi,Hy ) BUV, k) = 0. 



(43) 



This field equation is in general solved by a special 



function, namely the Whittaker function M\^^{x) 3]. We prefer to choose a particular 
linear combination of Whittaker functions, 



form of the confluent hypergeometric 



4' AniBHi 



+ /?^^^(fc)M(2) 



4' AttibH, 



■]k7] 



(44) 



with 

M(i)/(2) 



--e 4 



TL / nHir] 



r(|) 



k'^rf 



(45) 



because this tends to the asymptotic Bunch-Davies vacuum (i.e. the massless limit as r/ ^ 0), 



BUv.k) =aT{kW-H^P{kri)+^T{kW-Hf\k^) 

2 



r(2), 

2 



m. 



The Hankel functions are given by 

'r]^H^l\kr]) = 



(46) 



r]'^Hf\kr]) 

2 



Tik^ 



nk^ 



{—i — krj} e 



ikrj 



{i — kr]} e 



-ikrj 



With this choice of solution, the canonical condition, equation ()34|) is satisfied if 



W[BUri,k),B^:^{r,,k)] = {\P. 



rad|2 



, rad 1 2 
B 



(47) 
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The coefficients a and ^ in the solution are to be defined by a continuous matching of 
and its derivative at the time of transition from infiation to radiation era. For notational 
simphcity, we write M^^^ and M^^^ for the fundamental solutions in radiation era (|44j) . 
and -BjVij for the mode function 1)411) during infiation all evaluated at the infiation-radiation 
transition. This then leads to 



^ rad 



1 



M(i)M(2)' - M(2)M( 



Likewise we find 



Explicitly, that is 



i3. 



rad 
B 



Bt^m' -mB^:;\. 



^ rad 



rad 
B 



^ + o 

_1HI 
~2l^ 



Hj 



1 



k A;2 
2i 2^ 



me 



Thus we can see that for superhubble modes at the end of infiation, 

f^B 



jrad _ ^rad 



Hi' Hi 



2 A;2 



(48) 
(49) 

(50) 
(51) 

(52) 



This Bogoliubov mixing [IJ, of the modes is caused by the mismatch of the vacua during 
infiation and radiation era, and loosely speaking can be thought of as particle production 
at the infiation-radiation transition. 



C. Matter era 



During matter era the scale factor is given by 

The field equation for B^ during matter era reads 

6 , H] 



d'o + k''- -—-2 +<T^2(^ + ^e)' ) Bttiv, k) = 0. 



(53) 



(54) 
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Region 


Effective field equation 


Scale 


/ 




BLt{X) = 




II 




BLt{X) = 




III 




BLt{X) = 


X2 > ^'^^^'^ and X3 > 4^^=^^ 



TABLE I: The effective field equations for different values of X = kf] = 



A; (rye + f]). The expressions are written in terms of the parameter ^ 



As the field equation contains a term oc r^^, it cannot be solved analytically. However, for 
different values of krj the field equation effectively takes various forms, shown in tabelH] 
In region I the solution is given by 



2 



2 



(55) 



with again the coefficients given by a continuous matching at the radiation-matter equality. 



„ mat „ rad I a - i 



-ikrje _j_ ^I'^d 



In region // the solution of the field equation has to be solved numerically. In region ///, 



(56) 
(57) 



the field equation reduces to 



(a| + a;(r^)2)5l,(r^,/c) = 0. 



(58) 



If ^ < and ^ < u"^, which is the case for the demanded values of kfi in this region, the 
solution can be approximated by 

„mat _ /3mat 

obeying the condition 

1 



(59) 



Ui 



(60) 



The explicit expressions for a and (3 are rather complicated and for simplicity omitted. 
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V. THE SPECTRUM 



The power spectrum of the longitudinal mode, equation (|H7|). in momentum space reads 



In both radiation and matter era the subhubble spectrum grows as PsupcrH(^) oc k'^, 
which is an infrared save spectrum. In the ultraviolet region the spectrum converges to the 
Bunch-Davies vacuum, PsubH oc constant. In the intermediate region, the mass determines 
the height and blueness of the peak in the spectrum, hence in both radiation and matter era 
the mass term causes an enhanced spectrum with respect to the massless theory P|. The 
evolution of one mode, fixed k, in both radiation and matter era is given in figures ^ and |21 

Evolving in time, the damping of the modes in the density fluctuations evolves as "P oc 
on superhubble scales in both eras. For small enough mass, the spectrum remains relativistic 
in subhubble regions, that is, a damping of P oc a~^. However, as soon the mass term 
m%a{r])'^ becomes dominant in the field equation, the evolution of a mode becomes non- 
relativistic, V oc a~^, with a decaying oscillatory behaviour on top. A snapshot of the 
spectrum at fixed time r] is illustrated in figures El and El 

It is the aforementioned non-relativistic evolution, depending on the scale factor, which 
makes the theory viable for providing the energy density of dark matter. 



As the mass term in the field equations grows with the scale factor, at a certain mo- 
ment the energy density of the field switches from being dominated by relativistic modes to 
domination by non-relativistic modes, as previously explained. 

When this switch in domination occurs, the energy density of the field grows with respect 
to the total energy density of the universe during the radiation dominated era. 

If we assume that inflation can be described by a scalar inflationary model, such that 




(61) 



VI. ENERGY DENSITY 




3M2' 



(62) 
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FIG. 1: The evolution of a specific mode in momentum space, plotted for several val- 
ues of niBHi/k^. For low mass, a radiative spectrum is retrieved which redshifts oc 

a~^. The larger the mass, the earlier the evolution of a mode becomes massive as 
the mass exceeds the Hubble factor, where the redshift becomes oc (x r] x a~^. 

where Mp = (SttGat)^^ ^ 2.4 x 10^^ GeV is the reduced Planck mass, the energy density of 
the radiation field at the end of inflation equals 

PrAVend) = 3i/|M|. (63) 

The energy density of the nonsymmetric metric field is given by 

/r]k 
-^rB{v,k). (64) 

In order to calculate the observable energy density, we subtract the conformal part of it, 
which effectively means that we impose an ultraviolet cut-off in the power spectrum. The 
excess energy-density comes from the modes which crossed the Hubble radius during infla- 
tion. The smallest of these momenta is given by the physical momentum at the beginning 
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FIG. 2: The evolution of a specific mode in momentum space during the mat- 
ter era, plotted for several values of tubHj /k'^ . For low mass, a radiative spec- 
trum is retrieved which redshifts oc a~^. Again, the larger the mass, the earlier 
the evolution of a mode becomes massive, with a redshift of cx a^^ cc y. a~^. 




— me Hy 1)^ = f " 

" ' ' ' ' 

12 5 W 20 50 100 



FIG. 3: A snapshot of the spectrum during the radiation era at a fixed moment 
in time, plotted for several values of niBHrf'. On superhubble scales the spec- 
trum grows as /c^, whereas in the ultraviolet the spectrum becomes flat. Higher 
mass enhances the spectrum and blueshifts the peak, which lies at /c ~ \/7nBHi- 

of inflation, kphys\in = k/a{rjin) = Hf, while the largest momentum is given by the physical 
momentum at the end of inflation, kphys\end = k/a{r]end) = Hj, where r] = r]in [r] = r]end) 
denotes the conformal time at the beginning (end) of inflation. The latter momentum is 
also taken as the ultraviolet cut-off. 

As the variable over which is integrated is ki], the minimum and maximum value of the 
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FIG. 4: Snapshots of the spectrum during the matter era, plotted for several val- 
ues of ^ = iriBHrjl. Figures (a)-(d) represent mBHrjl = 10~^,10~^,1 and 10, re- 
spectively. The maximum value of the power spectrum still lies at A; ~ \/niBHi- 



domain of integration are given by ainairf) and aendO'-q respectively. Thus the energy density 
of the antisymmetric tensor field becomes 



Pb 



Hj 



1 



47r2a(r/)4 

where x = kr], = a{rjin) and a^nd 
that we may expand for small Om, 



cos 2x sin 2x 
log X H — I ; ci2a; 



2x 



(65) 



o-ijlend)- From equation (ji^ we have aend = 1, such 



Pb 



Hj 



log a(r7) - 1 + log2 + 7s + ]^al^a{r^f + O {a{r^)-^) + O (a^J 



(66) 



47r2a(r/)4 

Here 7^ denotes the Euler constant. Now we see that the energy density contains a subhubble 
term scaling as 0(77)"^, a term scaling as 0(77)"^ log 0(77) representing the continuous filling 
up of the density on horizon scales, and a small term aina{rj)~'^ for the superhubble energy 
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density. We assume the scale of inflation to be large enough to have ainairj) < 1, such that 
the expansion is indeed allowed. 

When the energy density is dominated by non-relativistic modes, the calculation of the 
energy density is less trivial. As the field is given in terms of M .^[—iniBHirf)^ the 

integral over all fourier modes too complicated to be evaluated exactly in an analytical 
fashion. However, we do know the asymptotic behaviour of the field in the infrared and the 
ultraviolet regions separately. 

In the infrared the leading order of the spectrum is oc [kr])'^. In the ultraviolet the 
spectrum becomes that of the massless theory. 

This information together with a closer look at the field equation, 

(^0^ + fc^ - ^ + mlHW^ BUv, k) = 0. (67) 

allows a simple approximation. Note that during the radiation era the term oc r^"^ represents 
a(?7)^if (r/)^, where H{ri) is the Hubble factor. 

The approximation is the following: each of the second, third and fourth term in equa- 
tion (jUTjl has its own region in which it defines the behaviour of the field, but for each value 
of k the mass eventually will dominate, leading to a scaling of a{ri)~^. Thus the energy 
density during the period of massive behaviour can be given in terms of the energy density 
during the period of massless behaviour. 

The superhorizon modes all together are the first ones to start scaling as non-relativistic 
matter, as soon as the Hubble constant shrinks about below the mass, trb > V2H{r]). That 
is the case when 0(77) ^ Hj/niB, such that 



PsuperH(A:,^) ^ nupSr(A:,r]) - a{v) ^ j^Vrj^^-{k,r^). (68) 

For sub-Hubble modes, the field starts scaling massive as soon as the wave vector shrinks 
below the mass, > k/aijj), when aijj) ^ k/ms, such that 

PsubH(A;, v) ^ VrXt-ik, ^^^^^Th ^^^(A:, v) (69) 

for H{ri) < < niB- The peak of the massive spectrum lies at the mode of which the 
physical momentum is equal to the mass when it crosses inside the Hubble radius, that is, 
k'^ = a{'qyH{viY = a{r))^m%, 

kpeak =mB a{riNR) 



-VHimB, (70) 
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where tinr denotes the moment in time of switching from relativistic to non-relativistic for 
the mode kp^ak- 

Now we can approximate the non-relativistic energy density by 



J TUB air]) 

^\b- 2a(r,)-log - 2, 



{k,r)) 



87r2a(r/)3Y Hj 



Hi V Hi 



+0 {ala{r)f, a{r))-', {^rriBHir))-')] . (71) 

In this expansion, B is the leading order coefficient, which can be determined by a numerical 
calculation of the energy density. In our approximation with a sudden relativistic-to-non- 

13 3 

relativistic transition, we find B ~ —24 -\- 2 + 2* log 2 + 247^ ~ 1.35. The third integral 
term only contributes if mBCi{T]) < Hi. Of course, when mBCi{r]) > Hi, the upper bound of 
the second integral becomes QeHi, and the logarithmic term in the energy density becomes 
nonexistent. This logarithmic term represents the contribution of superhubble but relativis- 
tic modes. This term damps quickly, and is neglected in the following. We assume that the 
scale of infiation is large enough to let aj„a(7y) -C 1 always, and we already made use of the 
fact that during the non-relativistic behaviour of the theory \/rnBHir] > 1. 

Taking account only for the leading order, we find for the NGT-to-radiation ratio when 
the field has become massive, 

Pb ( . fm^ r HfB 
Prad V Hi STimf, 

where nip — SnMp is the unreduced Planck mass. 

In order for the energy density to fully take account for the dark-matter energy density, 
it has to satisfy 

=1 at77 = 77e. (73) 



Prad 

Hence, 



(74) 
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This gives a rough approximation for the mass, 

Hi _ Heq 



TUB 



or 



m^^ ~ 7 X 10" 



Hi J 



Hi 



m. 



(75) 



(76) 



1013 GeV^ 

In the last hne we used for the Hubble parameter at radiation-matter equality, H^^ = (^eq + 
if'^Ho, Ho ~ 1.4 X 10-^2 Qgv, = 3230 ± 200 and the scale factor a^/ao = l/(^eq + 1)- 



VII. PRESSURE 



Similar to the energy density, the pressure of the theory can be calculated as the pressure 
is given by the spatial components of the diagonal of the stress-energy tensor. 



T,' = -3P. 



(77) 



1 ' I 1 2 

If isotropy is assumed, that is kik^ — ^Si^ \k\ for i — j, the expression for the 'pressure per 
mode' P becomes, 

2 



P{ri,k) 



2A7r^a{r]y 



doB\rj,k) + ^B\r),k) 
a{r]) 



-3k' 



B^{r),k) -mla{r)f B^{r),k) 



(78) 



The physical pressure can then be calculated by 

/rik 
-Piv,k). (79) 

When the field is dominated by relativistic modes, during radiation era this leads to. 



- log a(r7) + -log2 +37^-2 



1 2 . ^2 , sin2a(r/) 



+ 0{aia{vr,a{v)-') 



Covariant conservation of the stress-energy tensor implies 

a{r]) 



(80) 



(81) 
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FIG. 5: The 'pressure per mode' for non-relativistic modes during radiation era, for mBHi'rf = 
10^. Note that the X-axis is logarithmic, where the y-axis is Unear. For super hubble 
scales the pressure scales as —\p, whereas on subhubble scales is averages aroung zero. 

Comparing expressions and (jHUj) . we find that energy conservation is indeed obeyed 
order by order in rj. 

On superhubble scales, krj <^ 1, the power spectrum and the pressure per mode are 
dominated by the first terms in relations (j6T| and (|78|) . Thus one finds that for superhubble 
modes, both relativistic and non-relativistic, we have P{ri,k) ~ — |Pb(^5^), in accordance 
with the scaling of p oc 0(77)"^. 

During inflation and during the matter era similar results are found for the relativistic 
modes, satisfying energy conservation. 

This integration must be performed for the massive modes as well. However, for the 
same reasons as for the energy density, the pressure for the massive modes can only be 
either approximated or integrated numerically. We do not perform here a complete analysis 
of the pressure. Instead, we illustrate the pressure for one specific choice of parameters, and 
comment on how the pressure behaves for a general choice of parameters. 

In figure El the 'pressure per mode' is illustrated for a massive nonsymmetric field, with 
niBHirf = 10^. On superhubble scales, that is ki] <^ 1, the pressure per mode behaves 
as mentioned above as P{ri,k) = —^VB{v,k). On subhubble scales the pressure per mode 
becomes zero, conserving energy with a density scaling as p oc a{r])~^. 

During inflation and matter era the behaviour of the pressure is consistent with covariant 
energy conservation as well. The result is that we have zero pressure at the time of decoupling 
when the field is massive. This is in agreement with the requirements on cold dark matter. 
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VIII. DISCUSSION AND CONCLUSIONS 

We have considered the cosmology of a massive antisymmetric tensor field whose dynamics 
is given by the action (PJ. We follow the evolution of the vacuum fiuctuations generated 
during an infiationary epoch, through infiation, radiation and matter era. We find that the 
antisymmetric field with a mass, = 3 x 10^^ (10^'^ GeV/i//)^ eV, results in the right 
energy density today to account for the dark matter of the Universe. This then implies 
that below the corresponding length scale, m^^ ~ 0.1 /xm (////lO^^ GeV)^jthe strength of 
the gravitational coupling may change, as it has been argued in Refs. Q,^^. This scale is 
about two orders of magnitude below the present experimental bound, which is of the order 
10/xm [3]. Note that the mass scale (fTKjl depends strongly on the scale of infiation, such 
that, if the gravitational force law below the scale of 0.1 fim remains unchanged, would imply 
that, either the B field is nondynamical, or the infiationary scale is lower than 10^^ GeV. 
Note that, in contrast to gravitons, the antisymmetric tensor field begins scaling as non- 
relativistic matter during radiation era, making it potentially the most sensitive probe of 
the infiationary scale. 

The peak in the energy density power spectrum, generated as a consequence of a break- 
down of conformal invariance in radiation era, corresponds to a comoving momentum scale 
of the order, k ^ ^/mBHo/ {I + Z(,qY^^, where z^q — 3230 is the redshift at matter-radiation 
equality. Hence the most prominent matter density perturbations at a redshift z = 10, when 
structure formation beg ins, occur cit cl sccilc, ^p^ys 

~ 2 X 10^ km (corresponding to a wave- 
length, Aphys ~ 1 X 10* km, which is (co-)incidentally the Earth-Sun distance), which may 
boost early structure formation on these scales. We therefore expect that, when compared to 
other CDM models, this type of cold dark matter may induce an earlier structure formation. 
Even though the mass of the antisymmetric field is quite small (of the order of the heaviest 
neutrino mass, rrii, ~ 0.1 eV), the antisymmetric tensor field dark matter is neither hot nor 
warm. Indeed, since the antisymmetric tensor field does not couple to matter fields, it can- 
not thermalise and its spectrum remains primordial, and thus highly non-thermal. Hence, 
in spite of its small mass the field pressure remains small, such that structure formation on 
small scales does not get washed out. Indeed, P/pb — EkinI ijnBC^) ~ 10^^^ evaluated at 
kphys = {kpeak/a)\today, which is a tiny pressure. 

Although the pressure converges to zero as the field becomes dominated by non-relativistic 
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modes, the pressure per mode in Fourier space shows a characteristic oscillatory behaviour 
shown in figure El These oscillations may infiuence cosmological perturbations in a manner 
analogous to Sakharov oscillations [2^, which in turn may produce an observable imprint 
in the cosmic background photon fiuid. This question deserves further investigation. 
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